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COMBUSTION CONTROL SYSTEM AND METHOD FOR 
SPARK-IGNITION INTERNAL COMBUSTION ENGINE 

BACKGROUND OF THE INVENTION 
5 [0001] The present invention relates to a combustion control system and 
method for a spark-ignition internal combustion engine. 

[0002] A spark-ignition internal combustion engine is generally operated at 
a higher compression ratio for improvement in thermal efficiency, but such a 
performance improvement can be affected by engine knock. 
10 [0003] Hence, a knock control system that monitors the occurrence of 
engine knock with a knock sensor and performs feedback control on ignition 
timing has been proposed and already come into use so as to prevent engine knock 
while maximizing engine thermal efficiency. 

[0004] It is herein conceivable that the engine knock could be prevented 
15 more assuredly by not only carrying out a knock avoidance operation upon 

detection of the engine knock but also predicting knock occurrence according to a 
simulation model. 

[0005] Further, the ignition timing is commonly set using a prescribed 
ignition-timing map with reference to engine speed and power and then adjusted 

20 on detection of the engine knock. A multiplicity of actual vehicle experiments 
must be conducted to prepare the ignition-timing map, and the number of 
experiments required to prepare the ignition-timing map becomes increased with 
the recent application of variable mechanisms including a variable valve 
mechanism and an exhaust gas recirculation mechanism. This results in higher 

25 cost and greater facility for engine development. It can be however expected 
that the prediction of knock occurrence by a simulation model will allow an 
omission of many or at least some of the experiments for preparation of the 
ignition-timing map to thereby provide a great advantage in engine development. 
[0006] In view of the foregoing, Japanese Laid-Open Patent Publication No. 

30 7-332149 proposes a combustion control device for an internal combustion engine, 
capable of predicting the autoignition timing of an unburned air-fuel mixture, or 



equivalently, the occurrence of engine knock and controlling combustion so as to 
avoid the engine knock in accordance with the prediction. 

SUMMARY OF THE INVENTION 
[0007] When the engine knock occurs at a high intensity, there arise various 
problems such as a deterioration in engine durability and a unpleasant knocking 
sound. On the other hand, these problems do not come up when the engine 
knock occurs at a very low intensity. It is possible to operate the engine 
continuously under so-called trace knock (sporadic weak knock)" by setting the 
ignition timing. In other words, the engine combustion can be controlled 
properly by e.g. retarding the ignition timing in such a manner as to limit the 
knock intensity to a trace knock level. 

[0008] In the above-disclosed combustion control device, however, the 
nonoccurrence of autoignition is made a necessary condition for the knock 
avoidance without any consideration of the knock intensity. It is because there 
has not been enough analysis made on the correlation between the autoignition 
prediction and the actual knock phenomenon. As a result, the ignition timing 
becomes overretarded to cause an undesired drop in engine thermal efficiency. 
[0009] It is therefore an object of the present invention to provide a 
combustion control system and method for a spark-ignition internal combustion 
engine, by which the occurrence of engine knock can be prevented more properly 
based on the prediction of autoignition of an end gas, so as to maintain maximum 
engine thermal efficiency. 

[00010] According to a first aspect of the invention, there is provided a 
combustion control system for a spark-ignition internal combustion engine, the 
system being configured to: detect engine operating conditions; predict, based on 
the detected engine operating conditions, autoignition timing of an end gas and an 
amount of heat released due to autoignition of the end gas; and control 
combustion to establish such a relationship between the autoignition timing and 
the amount of heat released due to the autoignition as to give a knock intensity not 
higher than a specified intensity limit. 

[0001 1] According to a second aspect of the invention, there is provided a 
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combustion control system for a spark-ignition internal combustion engine, the 
system being configured to: detect engine operating conditions; predict, based on 
the detected engine operating conditions, an autoignition timing of an end gas and 
an amount of heat released due to autoignition of the end gas; calculate a knock 
5 intensity from the autoignition timing and the amount of heat released due to the 
autoignition; and control combustion in the engine in such a manner that the 
knock intensity is lower than or equal to a specified intensity limit. 
[0012] According to a third aspect of the invention, there is provided a 
combustion control method for a spark-ignition internal combustion engine, 

10 comprising: detecting engine operating conditions; predicting, based on the 

detected engine operating conditions, autoignition timing of an end gas and an 
amount of heat released due to autoignition of the end gas; and controlling 
combustion to establish such a relationship between the autoignition timing and 
the amount of heat released due to the autoignition as to give a knock intensity not 

15 higher than a specified intensity limit. 

[0013] According to a fourth embodiment of the invention, there is provided 
a combustion control method for a spark-ignition internal combustion engine, 
comprising: detecting engine operating conditions; predicting, based on the 
detected engine operating conditions, autoignition timing of an end gas and an 

20 amount of heat released due to autoignition of the end gas; calculating a knock 
intensity from the autoignition timing and the amount of heat released due to the 
autoignition; and controlling combustion in the engine in such a manner that the 
knock intensity is lower than or equal to a specified intensity limit. 
[0014] The other objects and features of the invention will also become 

25 understood from the following description. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0015] FIG 1 is a schematic view of a spark-ignition internal combustion 
engine according to one exemplary embodiment of the present invention. 
[0016] FIG 2 is a flowchart for controlling the operations of the 

30 spark-ignition engine according to one exemplary embodiment of the present 
invention. 



4 



[0017] FIGS. 3 A and 3B are diagrams of a two-zone combustion model. 
[0018] FIG. 4 A and 4B are graphs showing a relationship between knock 
intensity and autoignition timing and heat release amount under certain engine 
operating conditions. 
5 [0019] FIG. 5 is a graph showing variations in the autoignited mixture 
fraction with respect to the crank angle at the knock occurrence under different 
octane-number conditions. 

[0020] FIG. 6 is a graph showing variations in the autoignited mixture 
fraction with respect to the crank angle at the knock occurrence under different 
10 octane-number and compression-ratio conditions. 

[0021] FIG. 7 is a graph showing variations in the autoignited mixture 
faction with respect to the crank angle at the knock occurrence under different 
engine-speed conditions. 

[0022] FIGS. 8 A and 8B are diagrams of how to estimate a trace knock 
15 point on the prediction of autoignition timing and heat release amount according 
to one exemplary embodiment of the present invention. 

[0023] FIG. 9 is a graph showing the validity of a combustion simulation 
according to another exemplary embodiment of the present invention. 

DESCRIPTION OF THE EMBODIMENTS 
20 [0024] A first embodiment of the present invention will be now explained 
below with reference to FIGS. 1 to 8. 

[0025] FIG. 1 is a schematic view of a spark-ignition internal combustion 
engine according to the first embodiment. 

[0026] The engine has cylinder head 1, cylinder block 2, piston 3, 
25 combustion chamber 4 (defined by cylinder head 1, cylinder block 2 and piston 3), 
intake valve 5, valve timing control mechanism 5b, exhaust valve 6, intake port 7, 
exhaust port 8, fuel injection valve 9 and spark plug 10. Combustion chamber 4 
is supplied with fresh air through intake valve 5 and intake port 7. For air intake 
control, the open and close timing of intake valve 5 is adjusted with valve timing 
30 control mechanism 5b. Fuel injection valve 9 is disposed in intake port 7 to 
inject an appropriate amount of fuel and thereby form a flammable air-fuel 
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mixture in combustion chamber 4. The air-fuel mixture is ignited by spark plug 
1 0, and the burned product gas is discharged out of combustion chamber 4 
through exhaust valve 6 and exhaust port 8. The engine further includes engine 
control unit (ECU) 11 and various sensors, such as accelerator opening sensor 12, 
5 coolant temperature sensor 13, crank angle sensor 14, air flow meter 15, intake 
temperature sensor 16 and throttle sensor 17. These sensors 12 to 17 are 
connected to ECU 1 1 so that ECU 1 1 controls the operations of the engine 
integratedly in response to signals from sensors 12 to 17. 
[0027] FIG. 2 is a flowchart for controlling the operations of the engine 

10 according to the first embodiment. 

[0028] ECU 11 first receives the signal from sensors 12 to 17 to detect 
engine operating conditions including an accelerator opening, a coolant 
temperature, a crank angle (engine speed), an intake air amount, an intake air 
temperature and a throttle opening. By required signal processing and 

15 computing, ECU 1 1 estimates an effective compression ratio and valve timing for 
controlling valve timing control mechanism 5b, fuel injection valve 9 and spark 
plug 10. ECU 11 also estimates an air- fuel ratio, a residual gas rate and 
temperature and pressure in combustion chamber 4 at the close timing of intake 
valve 5. ECU 1 1 does a cycle simulation using the above estimated status values 

20 as initial parameters to predict autoignition timing of an end gas (unburned 

air-fuel mixture) and an amount of heat released due to autoignition of the end gas 
(hereinafter referred to as an "autoignition heat release amount"). The cycle 
simulation is done according to a two-zone combustion model and an autoignition 
model in the first embodiment. ECU 1 1 calculates a knock intensity from the 

25 predicted autoignition timing and heat release amount. ECU 11 further 

calculates, from the knock intensity, a trace knock point at which trace knock 
occurs and a minimum ignition advance for the best torque (MBT), i.e., ignition 
timing that presents optimal engine torque and specific fuel comsumption. 
Based on the trace knock point and MBT, ECU 1 1 performs combustion control in 

30 such a manner that the knock intensity becomes equal to or smaller than a trace 

knock level (an allowable limit of engine knock). ECU 1 1 gives feedback on the 
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trace knock point and MBT to set a heat release pattern in the cycle simulation. 
[0029] FIGS. 3 A and 3B diagrams of the two-zone combustion model, 
showing combustion chamber 4 when viewed from above. In FIGS. 3A and 3B, 
the flame front of a combustion flame is indicated by a heavy line, and the 
5 position of spark plug 1 0 is indicated by a cross. 

[0030] In the two-zone combustion model, it is assumed that the combustion 
flame propagates spherically from the position of spark plug 10 as shown in FIGS. 
3 A and 3 B, so as to divide combustion chamber 4 by the flame front into two 
zones: a burned zone in which the burned air-fuel mixture exists and an unburned 

10 zone in which the end gas exists. During the propagation of the combustion 

flame, the mass and enthalpy of such a thermodynamic system is transferred from 
the unburned zone into the burned zone. The unburned zone becomes 
compressed to a high-temperature and high-pressure state under the influence of 
the heat release and thermal expansion in the burned zone. 

15 [003 1] In order to predict the autoignition based on the temperature and 

pressure of the end gas, the autoignition model is then applied to the history of the 
end gas temperature and pressure data in the cycle simulation. The occurrence 
of engine knock is judged in the case where autoignition occurs in the unburned 
zone before the combustion flame has propagated throughout combustion 

20 chamber 4. The autoignition timing and heat release amount are derived upon 
judging the knock occurrence. 

[0032] The procedure of the cycle simulation by the two-zone combustion 
model and the autoignition model in the first embodiment will be next explained. 
In the following equations, the subscripts "Z>" and "w" are included to denote the 
25 burned zone and the unburned zone, respectively. 

[0033] The energy conservation in the burned and unburned zones per unit 
time are given by the equations (1) and (2), respectively: 

md{x b u b )+ pdV b -mdx b h u -dQ bwali = 0 (1) 

md({\-x b )u u ) + pdV u + mdx b h u -dQ uwall =0 (2) 

30 where m is the mass of the air-fuel mixture, x is the mass fraction, u is the specific 
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internal energy, i.e., the differential coefficient of internal energy,/? is the pressure, 
Vis the volume, dx is the change in mass fraction, h is the specific enthalpy and Q 
is the amount of heat transfer. The third term mdxbh u of the equations (1) and (2) 
represents the amount of enthalpy transferred from the unbumed zone to the 
5 burned zone via the flame front, and the fourth terms dQb^aii and dQ UyWa ii of the 
equations (1) and (2) represent the amounts of heat transferred from the wall of 
combustion chamber 4 into the burned zone and the unburned zone, respectively. 
[0034] The equation (3) holds on the assumption that combustion chamber 4 
is a closed system where the change in volume per unit time is preserved 
10 according to the two-zone combustion model: 

dV b +dV=dV (3). 
The total volume Kof combustion chamber 4 may be given as a function of the 
crank angle ft 

[0035] The following characteristic equations (4) and (5) hold for the 

15 burned zone and the unburned zone, respectively: 

^_ + dV^^^ + dR L+ dT^ 
P V b x b R b T b 

p V u l-x b R u T u 
where R is the gas constant and Tis the gas temperature. 
[0036] Further, the mass fraction x b of the burned gas is given by the 
20 following Wiebe function (6) in the first embodiment: 

v m+1 



x b =1- exp 



-a 



0-Oq 
A0 



(6) 



where 0 0 is the crank angle of spark ignition, AO is the crank angle range during 
the combustion, and a and m are model constants. Alternatively, the mass 
fraction x b of the burned gas may be retrieved from an experimentally 
25 predetermined heat release pattern table. 

[0037] The changes of states in the burned and unburned zones are 
estimated by solving the simultaneous equations (1) to (6) for the crank angle 0. 
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To solve the equations (1) to (6), the initial parameters are set to the above 
estimated status values. Also, both the fourth terms dQ DtWa u and dQ UM , a ti of the 
equations (1) and (2) are taken as zero on the assumption that the cylinder gas is 
thermally insulated from the wall of combustion chamber 4. 
5 [0038] The occurrence of engine knock is judged depending on whether or 
not the end gas is autoignited before the combustion flame has propagated 
throughout combustion chamber 4 according to the autoignition model as 
mentioned above. Then, the ignition delay x of the end gas is calculated to 
predict the autoignition timing and heat released amount. 
10 [0039] In the first embodiment, the ignition delay x of the end gas under a 
certain high-temperature and high-pressure state is given by the following 
equation (7): 

r = Ap- expj^|j (7) 

where A, B and n are model constants. For prediction accuracy, the equation (7) 
15 needs to be applied where the end gas has already been compressed to a certain 
high-temperature high-pressure state due to increase in the cylinder pressure. It 
is thus most appropriate to apply the equation (7) to a crank angle Q pmax at which 
the cylinder pressure becomes maximized in the course of solving the equations 
(1) to (6) for the crank angle 0. 
20 [0040] The end gas ignition delay x at the crank angle Q pmax is converted 
from the unit time to the unit angle and added to the crank angle Q pmax , thereby 
giving a predicted value Q knock of the autoignition timing. Further, the predicted 
autoignition crank angle Qknock is substituted into the equation (6) to determine the 
current burned mixture fraction, i.e., the current heat release amount. The 
25 current burned mixture fraction is subtracted from 1 (one), thereby giving the 

unburned mixture fraction at the crank angle Qknock, i.e., a predicted value X a uto-ign 
of the autoignition heat release amount. 

[0041] Herein, there is a specific relationship between the knock intensity 
and the autoignition timing and heat release amount as discussed below. 
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[0042] FIG 4A is a graph showing variations in the normalized burned mass 
rate Z,„ (deg.* 1 ) with respect to the crank angle 0 (deg. ABDC), and FIG. 4B is a 
graph showing variations in the autoignited mixture fraction with respect to the 
crank angle at the knock occurrence Qknock (deg. ABDC). It is noted that FIGS. 
5 4A and 4B are quoted from SAE Technical Paper Series, Yasuo Takagi et al., "An 
Analytical Study on Knocking Heat Release and its Control in a Spark Ignition 
Engine," No. 880196, pp. 1-10, 1988. 

[0043] As is generally known, the autoignition heat release amount 
increases with the spark advance toward top dead center (TDC), resulting in a 

10 higher knock intensity, when the engine operating conditions other than the 

ignition timing are maintained. Such a relationship between the ignition timing, 
the autoignition heat release amount and the knock intensity under certain engine 
operating conditions is clearly proved in FIGS. 4A and 4B, but it has been unclear 
whether the same relationship can be universally established even when not only 

15 the ignition timing but also the other engine operating condition or conditions, 
such as research octane number and compression ratio, are changed. 
[0044] Accordingly, the present inventor has conducted the analysis of 
experiments on engine knock under various conditions to examine and generalize 
the relationship between the knock intensity and the autoignition timing and heat 

20 release amount. The analysis gives the following findings. 

[0045] FIG 5 is a graph showing variations in the autoignited gas mixture 
fraction (%) with respect to the crank angle at the knock occurrence Qknock (deg. 
ATDC) under different octane-number conditions. 

[0046] When the research octane number RON is 90, the autoignited gas 
25 mixture fraction varies with the crank angle Qknock in qualitatively and 

quantitatively the same way as in FIG 4B. On the other hand, the autoignited 
gas mixture fraction stays substantially the same with respect to the crank angle 
Qknock when the research octane number RON is 100. The crank angle Qknock 
corresponds to the autoignition timing, while the autoignited gas mixture fraction 
30 corresponds to the autoignition heat release amount. It is thus concluded that a 
higher octane number results in a larger shift of the autoignition timing toward 
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TDC and a smaller autoignition heat release amount, whereby the engine knock 
becomes more unlikely to occur. 

[0047] FIG 6 is a graph showing variations in the autoignited gas mixture 
fraction (%) relative to the crank angle at the knock occurrence Qknock (deg. 
5 ATDC) under different octane-number and compression-ratio conditions, obtained 
by averaging of the data of FIG. 5 on the cycle number. In FIG. 6, the data at a 
higher octane number and a lower compression ratio is plotted on the lower left 
side, and the data at a lower octane number and a higher compression ratio is 
plotted on the upper right side. 

10 [0048] At the same knock intensity, a larger advance of the autoignition 

timing results in a smaller autoignition heat release amount. In other words, the 
knock intensity increases with the autoignition heat release amount. As the 
autoignition timing is advanced, the knock intensity becomes relatively high even 
when the autoignition heat release amount is small. It is important to consider 

15 such a relationship between the knock intensity and the autoignition timing and 

heat release amount in the case where the knock occurrence is judged based on the 
prediction of end gas autoignition. 

[0049] FIG 7 is a graph showing variations in the autoignited gas mixture 
fraction (%) relative to the crank angle at the knock occurrence Qknock (deg. ATDC) 

20 under different engine-speed conditions. The data of FIG 7 at each engine speed 
represents trace knock points (at which the knock sensor signal K.L is 4). 
[0050] The correlation between the autoignition timing and heat release rate 
at each engine speed can be well approximated to an exponential curve. Even 
when the autoignition heat release amount decreases with increase in the engine 

25 speed, the engine knock occurs at the same intensity. In other words, the knock 
intensity increases with the engine speed. 

[005 1] Incidentally, a correlation curve of trace knock points at an engine 
speed between 1200 rpm and at 2400 rpm can be obtained by interpolation of the 
exponential correlation curves at 1200 rpm and at 2400 rpm. It is also possible 
30 to obtain a correlation curve of trace knock points at different knock intensity. 
[0052] It has been thus experimentally proved that the knock intensity 



11 



increases with the autoignition heat release rate, the advance of the autoignition 
timing and the engine speed. Based on such a universal relationship, the knock 
intensity is calculated from the predicted autoignition timing Q knock and heat 
release amount Xauto-ign in the cycle simulation. 
5 [0053] In the actual spark ignition engine, there is no need to perform a 
knock avoidance operation when the knock limit is on an advance side of MBT 
ignition timing, and it is most appropriate to set the spark ignition timing to the 
MBT ignition timing in view of the thermal efficiency. In the earlier technology, 
the MBT ignition timing is generally retrieved from a prescribed ignition-timing 
10 map with respect to the engine speed and load. However, the MBT ignition 

timing varies with the engine operating conditions, and there is a limit in setting 
the ignition timing using the ignition-timing map. 

[0054] It is thus reasonable in the first embodiment to simulate both of MBT 
and the trace knock point as engine torque can be computed in the cycle 

15 simulation concurrently with the knock intensity calculation. Namely, MBT and 
the trace knock point are determined by estimating the knock intensity and engine 
torque with respect to the heat release pattern based on the Wiebe function (6) and 
repeating, while shifting the heat release pattern in timing, the estimation of the 
knock intensity and engine torque over several or dozens of times. 

20 [0055] FIGS. 8A and 8B are schematic illustrations of how to the trace 

knock point. A curve a of FIG. 8B is a plot of the predicted autoignition timing 
Qknock and heat release amount X auto , ign obtained by the cycle simulation, whereas a 
curve P of FIG. 8B is a correlation curve: y = 0.826exp(0. Wlx) depicting the trace 
knock level. 

25 [0056] The trace knock point is defined as a point of intersection of the 

curves a and (3 of FIG. 8B, and then, converted into ignition timing presents trace 
knock (referred to trace knock ignition timing) as shown in FIG. 8 A. 
[0057] MBT is determined by simply finding through the cycle simulation 
an ignition timing point for the maximum engine torque. 

30 [0058] In order for the engine to maximize engine thermal efficiency while 
avoiding the occurrence of heavy engine knock that can affect engine durability, 
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the spark ignition timing is set to either one of the MBT ignition timing and the 
trace knock ignition timing located on a retard side. 

[0059] As described above, ECU 1 1 predicts the autoignition timing and 
heat release amount based on the engine operating conditions, and then, controls 
5 the combustion to establish such a relationship between the autoignition timing 
and the autoignition heat release amount that the knock intensity is lower than or 
equal to a specified intensity limit corresponding to the trace knock level. 
Further, ECU 1 1 calculates the knock intensity from the predicted autoignition 
timing and heat release amount so as to limit the knock intensity to the trace 
10 knock level assuredly. By adjusting the ignition timing in view of the knock 

intensity, it becomes possible to operate the engine without causing an undesired 
deterioration in engine thermal efficiency due to the overretardation of the ignition 
timing and thereby possible to maintain high engine thermal efficiency and 
durability. 

15 [0060] In the first embodiment, the autoignition timing and heat release 
amount are predicted on the ignition delay x of the end gas. In addition, the 
knock intensity is calculated in such a manner that the knock intensity increases 
with the autoignition heat release amount, the advance of the autoignition timing 
and the engine speed. This makes it possible to simplify the prediction of the 

20 autoignition timing and heat release amount and facilitate the combustion control. 
As the combustion control is performed by adjusting the spark ignition timing in 
the first embodiment, the time between the knock occurrence prediction and the 
actual combustion control can be reduced to avoid the engine knock more 
assuredly. 

25 [0061] Next, a second embodiment of the present invention will be 

explained below. The second embodiment is similar to the first embodiment, 
except that the cycle simulation is done in consideration of the temporal 
characteristics of the end gas ignition delay x as well as the heat transfer from the 
wall of combustion chamber 4 for high prediction accuracy. 

30 [0062] In the second embodiment, the following Livengood-Wu integral (8) 
is applied to the end gas ignition delay x: 



I. 
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The Livengood-Wu integral (8) is based on the concept that the end gas ignition 
delay x varies with the passage of time. 

[0063] According to the Livengood-Wo integral (8), the occurrence of end 
5 gas autoignition is judged at the time the integral of the inverse of the end gas 

ignition delay x with respect to the time t becomes equal to 1 (one). The end gas 
ignition delay x of the Livengood-Wu integral (8) may be expressed as a function 
of the temperature T and pressure p as in the equation (7). It is alternatively 
possible to store a preset ignition-delay map in ECU 1 1 and retrieve the end gas 
10 ignition delay x from the ignition-delay map with respect the temperature T and 
the pressure p. 

[0064] With the temporal characteristics of the end gas ignition delay x 
reflected in the cycle simulation as described above, it becomes possible to insure 
high prediction accuracy even when there arises relatively large temperature and 

15 pressure changes during the end gas ignition delay x. 

[0065] In order to reflect the heat transfer from the wall of combustion 
chamber 4 in the cycle simulation and thereby insure higher prediction accuracy, a 
known heat transfer model for the internal combustion engine is further applied. 
In the second embodiment, the so-called "Woschni's heat transfer model" is 

20 adopted. 

[0066] According to the Woschni's heat transfer model, the heat transfer 
from the wall of combustion chamber 4 is given by the following equation (9): 

_ dQwatl ~ A K ft - Tymll ) (9) 

where A is the heat transfer area, h c is the heat transfer coefficient and T wa u is the 
25 temperature of the combustion chamber wall. 

[0067] The heat transfer coefficient h c is given by the following equations 
(10) and (1 1) according to the Woschni's heat transfer model: 

h c = CB n '- l p m w m T 075 - l - 62m (10) 
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M,=c x T p +cJ*Lip- Pm ) (ii) 

PrK 

where C, G, C2, and m are the model constants, B is the bore diameter, p is the 

cylinder pressure, T is the average in-cylinder gas temperature, S is the 

average piston speed, Vd is the displacement volume, p m is the cylinder pressure at 
5 the motoring and the subscript "r" is included to denote the reference point set to 
e.g. the close timing of intake valve 5. The model constants Ci and C2 are 
changed in accordance to the phase of cycles of compression, combustion and 
expansion. 

[0068] The Woschni's heat transfer model is originally designed to address 

10 average heat transfer throughout combustion chamber 4. To combine the 

Woschni's heat transfer model with the two-zone combustion model, it is thus 
necessary to allocate the amount of heat transferred from the wall of combustion 
chamber 4 to the burned and unburned zones in some way. In the second 
embodiment, it is assumed that the end gas is in a thermally insulated state during 

15 the heat release process. Based on such an assumption, the fourth term dQ UtWa u 
of the equation (2) is taken as 0 (zero), and the fourth term dQb Ava ii of the equation 
(1) is taken as the heat transfer amount of the equation (9). 
[0069] The progression of the end gas temperature T u and pressure p u over 
time are thus estimated by solving the equations (1) to (6) and (8) to (1 1) 

20 simultaneously for the crank angle 0. In the course of solving the simultaneous 
equation (1) to (6) and (8) to (11), the occurrence of end gas autoignition is judged 
by calculating the end gas ignition delay x from the estimated values, and then, 
integrating the inverse of the calculated end gas ignition delay t with respect to 
time The predicted value Qknock of the autoignition timing is given as a crank 

25 angle at which the integral becomes equal to 1 (one) as in the equation (8). The 
predicted value X au to-ign of the autoignition heat release amount is given by 
substituting the predicted crank angle & knock into the equation (6) and then 
subtracting the thus-obtained current heat release amount from 1 (one). 
[0070] The knock intensity, the MBT ignition timing and the trace knock 
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ignition timing are calculated from the predicted autoignition timing Q knock and 
heat release amount X Q uto-ign in the same manner as in the first embodiment. 
[0071] As described above, the autoignition timing and heat release amount 
can be predicted accurately in view of the temporal characteristics of the end gas 
5 ignition delay x and the heat transfer from the wall of combustion chamber 4. 
This makes it possible to limit the knock intensity to the trace knock level more 
assuredly so as not to cause an undesired deterioration in engine thermal 
efficiency due to the overretardation of the ignition timing and thereby possible to 
maintain high engine thermal efficiency and durability. 

10 [0072] Finally, a third embodiment of the present invention will be 

explained. The third embodiment is similar to the first and second embodiments, 
except that the occurrence of end gas autoignition is predicted according to 
another autoignition model called "elementary reaction model". 
[0073] The elementary reaction model is based on the chemical kinetics in 

15 which an overall reaction is made up of a series of several tens to hundreds of 
elementary reactions so as to be conformable to actual physical phenomena 
(including molecule-to-molecule collision and reaction). The influence of not 
only the temperature and pressure but also the gas mixture composition and the 
fuel properties over the autoignition can be simulated by such an elementary 

20 reaction model. With the elementary reaction model applied to the end gas in the 
unburned zone, it becomes possible to predict the knock occurrence more 
accurately in response to the changes in valve timing, compression ratio, 
equivalence ratio, residual gas rate and octane number etc. 
[0074] The primitive equations (1) to (6) and the Woschni's heat transfer 

25 model equations (9) to (1 1) are similarly used in the third embodiment. 

[0075] It is now assumed that the gas mixture of combustion chamber 4 
contains a mass fraction^, of the chemical species i and that the mass fraction^, 
of the chemical species i varies during the reaction. 

[0076] The change in mass fraction y t of the chemical species i during the 
30 reaction is thus given by the following equation (12) according to the elementary 
reaction model: 
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dy, u =^o>, u (12) 

Pu 

where cb ,-,„ is the reaction rate of the species /, M t is the molecular weight of the 
species i and p w is the density. With the application of the Arrhenius equation to 
the reaction rate cb , iU , the reaction rate cb itU is expressed as a function of/?, T u 
5 andj/*,,,. Further, the reaction rate cb i>u can be be easily computed by 

CHEMKIN (a computer program commonly known for calculation of known 
homogeneous reaction kinetics). 

[0077] The change in internal energy du lt during the reaction is given by the 
following equation (13) according to the elementary reaction model: 

(I3) 

Pu 

[0078] Also, the change in gas constant R u during the reaction is given by 
the following equations (14) and (15) according to the elementary reaction model: 

<*K=-^frdM u (14) 
„ M„d>, „M, 

dM u = Z " ~ ' (15) 

/ Pu 

15 where R urU is the general gas constant and M u is the average molecular weight. 
[0079] Although not discussed in the specification, it is also possible to 
apply the elementary reaction model to the burned zone so as to simulate not only 
engine knock, i.e., the end gas autoignition in the unburned zone but also the 
thermal dissociation and NO generation in the burned zone. 

20 [0080] The progression of the end gas temperature T u and pressure p u and 
the mass fraction M t of the chemical species i over time is thus estimated by 
solving the equations (1) to (6) and (9) to (15) simultaneously for the crank angle 
0. As the cycle simulation on the elementary reaction model deals with a 
complex system that involves sudden state changes, a high-precision, stable solver 

25 needs to be used to solve the equations (1) to (6) and (9) to (15). 

[0081] The autoignition is regarded as a part of the continuous chemical 
reaction in the elementary reaction model. It is thus necessary to set a criterion 
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forjudging the occurrence of the autoignition. For example, the occurrence of 
the autoignition may judged when the end gas temperature T u becomes equal to or 
larger than a certain value. 

[0082] The predicted value Qknock of the autoignition timing is given upon 
5 judging the occurrence of autoignition based on the set criterion. Then, the 
predicted value X au to-ign of the autoignition heat release amount is given by 
substituting the predicted crank angle Qknock into the equation (6) and subtracting 
the thus-obtained current heat release amount from 1 (one). 
[0083] The knock intensity, the MBT ignition timing and the trace knock 
10 ignition timing are calculated from the predicted autoignition timing Qknock and 
heat release amount Xauto-ign in the same manner as in the first or second 
embodiment. 

[0084] FIG. 9 is a graph showing how the spark advance AADV varies with 
the compression ratio s in the actual vehicle experiment and in the cycle 
15 simulation of the third embodiment. As shown in FIG. 9, the prediction results 
well agree with the experimental results qualitatively and quantitatively at each 
engine speed. The validity of the prediction in the third embodiment is thus 
clearly proved in FIG. 9. 

[0085] As described above, the autoignition timing and heat release amount 
20 can be predicted more accurately in response to the changes in valve timing, 
compression ratio, equivalence ratio, residual gas rate and octane number etc. 
because the influence of not only the temperature and pressure but also the gas 
mixture composition and the fuel properties over the end gas autoignition is 
reflected in the cycle simulation by the elementary reaction model. This makes 
25 it possible to limit the knock intensity to the trace knock level more assuredly so 
as not to cause an undesired deterioration in engine thermal efficiency due to the 
overretardation of the ignition timing and thereby possible to maintain high 
engine thermal efficiency and durability. 

[0086] The entire contents of Japanese Patent Application No. 2003-126869 
30 (filed on May 2, 2003) are herein incorporated by reference. 

[0087] Although the present invention has been described with reference to 
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specific embodiments of the invention, the invention is not limited to the 
above-described embodiments. Various modification and variation of the 
embodiments described above will occur to those skilled in the art in light of the 
above teaching. The scope of the invention is defined with reference to the 
5 following claims. 



